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A fermion ground state energy functional is set up in terms of particle density, relative pair density, 
and kinetic energy tensor density. It satisfies a minimum principle if constrained by a complete set 
of compatibility conditions. A partial set, which thereby results in a lower bound energy under 
minimization, is obtained from the solution of model systems, as well as a small number of exact 
sum rules. Prototypical application is made to several one-dimensional spinless non-interacting 
models. The effectiveness of "atomic" constraints on model "molecules" is observed, as well as the 
structure of systems with only finitely many bound states. 
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I. INTRODUCTION 

The ground state properties of many- fermion systems, 
even with the implied restrictions on temperature do- 
main, stationarity, homogeneity of particle type, have a 
physical importance that can hardly be overestimated. 
This has spawned a large variety of approximative, and 
in principle exact analytical and numerical solution tech- 
niques. There are also, in practice, conceptual restric- 
tions on such techniques. They should not supply appre- 
ciably more information (which costs more as well) than 
one can conceivably use, but the intermediate constructs 
employed should be sufficiently nuanced that qualita- 
tive distinction between distinct systems can emerge, and 
"physical intuition" utilized and developed along the way. 

One "traditional" basically analytical direction that 
has been, and is being, pursued intensively is based upon 
the Rayleigh - Ritz variational principle. It is extreme 
in both aspects alluded to above: the intermediate con- 
struct used is the full A^-body trial wave function -0Ar , and 
the exactly bounded information obtained is the ground 
state energy £'o alone: 



-Bo < {tpN\HN 1-0^) 



(1) 



here, Hn denotes the system Hamiltonian. To be sure, 
bounding principles for expectations are also available 
but are much more conservative, the ground state restric- 
tion is trivially weakened to include ground state condi- 
tioned by compatible constants of the motion, and highly 
accurate ground state energy engenders confidence that 
expectations using the computed wave function will be 
highly reliable. And of course, the generality of this ap- 
proach makes application to multi-species systems quite 
direct, but also very computation-intensive. 

As it has become necessary to deal with increasingly 
complex structures, the electronic component of a Born - 
Oppenheimer macromolecule being an outstanding ex- 
ample, density functional techniques have emerged as 
very powerful semi-analytic tools. They use as intermedi- 



ate constructs, and produce as output, the same quantity, 
the fermion (electron) density 



N 



Pn{x) = {y25{xj -x) 

= {^Hx)i,{x)) , 



(2) 



in first and second quantized form respectively, x denot- 
ing all degrees of freedom of a particle. More elaborate 
versions work instead with the reduced one - body den- 
sity matrix 



{x\-in\x') = {^1^\x)^{x')). 



(3) 



and attend to the spin population as well, but all depend 
upon intelligent semi-empirical information as to the de- 
pendence of the energy on this quantity, corresponding to 
a bounding principle of not much more than pictorial sig- 
nificance. In practice, even primitive local density func- 
tionals often give tolerable results, and density gradient 
extensions jl| can, in professional hands, be remarkably 
effective, as well as be generalizable to to non-ground 
states and dynamical information. 

An elusive goal for many years has been the estab- 
lishment of practical analytical techniques which work 
through not much more than the needed output informa- 
tion and yet are also exact bounding principles at some 
level. The modifier "practical" is where the difficulty 
lies. The version that has been pursued even longer than 
density functional techniques is that involving only the 
ground state 2-body reduced density matrix 



{x^,x2^TN\x'^,x'^) = {^\x2W{xlmx'lmA)) (4) 



of a pair-interacting system of the form 
Hn = hi + ft,2, 

N 

i=l 

/l2 = ^ u{xi,Xj). 
l<i<j<N 

Then one has, in obvious notation, 



Eo= min Tr(X7Ar) + Tr(urjv), 

7jv,rjv 

where K — T -\- v . 



(5) 



(6) 



But the reduction of an TV-body system to the de- 
tailed properties of an effective 2-body system involves a 
tremendous reduction of information, which must there- 
fore be supplied indirectly [2|. In the version referred 
to, these take the form of a small selection of known 
sum rules, and a large implicit, but only partially known, 
selection of inequalities that the 2 body density matrix 
must satisfy. It is only with the advent of modern com- 
putational tools - hardware and algorithms - that this 
approach has become feasible, but still primarily in the 
area of small systems and coarsely discretized function 
space. 

In the present paper, we study an offshoot of this last 
activity, initially restricted to the very large domain of 
systems that can be modeled by having Newtonian ki- 
netic energy 



T{p)=py2r 



(7) 



as the only momentum - dependent contribution. The in- 
termediate quantities we deal with are various one-point 
densities, particle density jSJ as well as suitably defined 
kinetic energy tensor density and pair interacting density 
(The pair distribution integrated over center of mass). 
These are minimally informative, but sufficient to repre- 
sent major quantities of physical interest - including the 
energy, thereby allowing an energy bounding principle to 
be formulated. The reduction in information is even more 
extreme than in the pair-density matrix formulation, and 
so even more subsidiary restrictions must be imposed. 
But we will see that they can be selected by attention 
to the underlying physics, and then are rewardingly ef- 
fective. This means that we will take great advantage 
of simple physical systems that share physical charac- 
teristics with the system at hand. An important tool 
will stem from the observation that any solvable model - 
each of which is associated with an inequality restriction 
- expands to a whole class of models under coordinate 
transformations. 

In our current entree to this approach, we will confine 
our attention to very primitive "toy" models, selected 
to probe the effectiveness of the technique without 
undue complexity. Thus, starting in Sec II, we deal only 
with one-dimensional non- interacting spinless fermions. 



establishing the basic sum rules and the model-generated 
form of inequalities. As will be reported in work now 
in progress, extension to 3 dimensions and spin is quite 
direct, if a bit more complicted, whereas the explicit 
inclusion of interaction requires a simple but non-trivial 
expansion of the technique (not even needed if only 
the construction of the Kohn-Sham pseudo-one-body 
density matrix |3| is at issue). In particular, the required 
amassing of model reference systems is more intricate, 
and will be reported on in a later publication; see 
however Section V. In Sec III, we apply the technique 
to a few very elementary examples, and indicate in Sec 
IV how the strict minimization can be relaxed. Sec V 
tests to what extent toy "molecules" can be solved via 
knowledge of their component "atoms" , and in Sec VI, 
we study the much more demanding situation in which 
only a finite number of non-interacting bound state 
exists. 



II. IMPLICIT DENSITY FUNCTIONAL 
INEQUALITY TAKEN AS CONSTRAINT 

Since one knows, e.g. from the work of Hohenberg 
and Kohn 4], that one can use the A^-electron density 
function pn{x) as the only variable when interaction and 
external potentials u and v are fixed, let us rewrite Eq. ^ 
as 



Eo = min <^ T[pn] 

pn y 



PN{x)v{x)dx + f/intiPw] f ; (8) 



where T[pn] is the kinetic energy part and U-y^i^[pM] 
the interacting part of the energy expectation. The in- 
principle separation into T[piq] + U-y^^[pN] is a conse- 
quence of the fact that pN determines v and hence the 
full wave function -^^r as well; T[p7v] is not identical to 
the adiabatically interactionless kinetic energy relevant 
to Kohn-Sham. 

Therefore, the ground state can be obtained by apply- 
ing a variational principle to the expected value of energy 
with respect to the density function pn{x), provided that 
the full functional of density is known. In most cases, the 
exact form of the functional is impossible to write down 
explicitly. And as we will emphasize in this paper, it's 
also unnecessary to do so since the functional is implied 
by sufficient many equalities or inequalities, several of 
which we apply as constraints on the minimization. 

Consider a one dimensional non-interacting A'^-fcrmion 
reference system on coordinate space {A}, the ground 
state energy E^ of which is known and must satisfy 



{iPnIHn \iPn) > En, 



(9) 



where Hn — ^i=i{^Pi + ^{^i)), and V'a' is any anti- 
symmetric A'^-body wave function. An enormous conve- 
nience is that each solved reference generates a whole 



class of useful references, as follows: A continuous trans- 
formation of coordinates Xi = f{xj) can always be ex- 
tended to a unitary tranformationQ by setting 



B. 



g{xi)pi +Pig{xi) 



(10) 



where g{x) = \/J{x) and J{x) = f'{x) is the trans- 
formation Jacobian. Hence the commutator relation 
[xi,pj\ — ihSij persists, and a new set of canonical coor- 
dinates {xi,pi} is obtained. Define the symmetric kinetic 
energy density operator (which becomes a tensor density 
in higher dimensional space) as 

1 ^ 
f{X) = -Y, (KX - Xi)P,P, + 2P,5{X - Xi)Pi 



P^P^S{X - X,] 



(11) 



Ea. (|ll|l is not a unique representation of kinetic energy 
density, since we can add any spatial divergence to it 
and yield the same total energy. For example, one has 
equivalent symmetrized up-section 



Tu{X) 



1 ^ 

i=l 



X,),P,P, 



(12) 



or mid-section 



tn{X) 



1 ^ 

-Y,mx-Xi)p,. 



(13) 



i=l 



Both lead to the same total kinetic energy but the ki- 
netic energy density is different. For a non- interacting 
fcrmion system with harmonic oscillator external poten- 
tial well, the three types of kinetic energy densities men- 
tioned above are shown in Figure ^ Compared to the 
other two, the shell details of the fully symmetrized ki- 
netic energy density are diminished. The advantage of 
picking up the symmetric one is that it gives a relative 
simpler form of the coordinate transform for kinetic en- 
ergy. 

Applying the symmetrized kinetic energy density def- 
inition Hll|l . we have for the energy expectation for the 
previous Hamiltonian 

{iPN\HN\i^N)^ fT{X)dX+ f pv{X)V{X)dX 



= I g'^{x)t{x)dx + I py{x)v{x)dx (14) 
+ o / Pv{x)g''^{x)dx, 



where T and pv are kinetic energy and electron density 
for the exactly solved reference system, t{x) is the trans- 
formed kinetic energy density obtained with the symmet- 
ric kinetic energy density operator, defined the same way 



as T{x): 
i{x) 



1 ^ 

- ^ {S{x - Xi)piPi + 2piS{x - Xi)pi 



+ PtPiS{x - Xi)) 



(15) 



Moreover, v{x) = V{X) and p„ is transformed from the 
original density function as Pvix) = pviX)J{x). 

Therefore, for any realizable, i.e. "iV-representable" 
combination of kinetic energy density t(x) and density 
p{x), we must have 



g^ix)t{x)dx+ p{x)V{f{x)) 



(16) 



p{x)g'^{x)dx>EN[V], 



where f'{x) = J{x) — \/ g{x). Since the transform J{x) 
is arbitrary, we have an inequality to be satisfied by any 
A'^-representable combination oit{x) and p{x): 



(17) 



With V{X) and E]si\V] given by the reference system, 
we have a well defined constraint for the functional re- 
lation between t(x) and p{x). A lower energy bound is 



min / g {x)t{x)dx+ / p{x)V {f (x)) dx 

+1 f Pix)g'\x)dx > En[V] 



•symmetric 

-mid-section 

- up-section 




-4 



-2 




X 



FIG. 1: (Color online) Three types of kinetic energy density for 
the ground state of harmonic system with 6 non-interacting 
fermions are shown here. The thickest dark line represents 
the symmetrized kinetic energy density, while the other two 
are the up- and mid-section of the symmetric kinetic energy 
density (Eq.(|nj)- 



obtained by carrying out the minimum with only the con- 
straint satisfied, because we certainly have not included 
all constraints needed to guarantee A^-representability. 

For a given v-representable (belonging to the potential 
v) density function Pv(x), according to inequality (|16|l . 
the coordinate transform leads to 



F\V] = f gv\x)t{x)dx + ^ J p{x)gv'\x)dx ~ En[V] 

> — / pi,{x)v{x)dx, 

(18) 

where gv{x) = V (y~^ {v (x))) /v'{x). As this inequal- 
ity is ubiquitous for any V{X), we must have 



minF[y] -|- / py{x)v{x)dx > 0. 



(19) 



According to Ea. (|18|l . it's obvious that the minimum for 
F[V] will be achieved if g{x) = 1 where the inequality 
Eq. (|19|l becomes an equality. Approaching its minimum, 
we have tttW = 0, which gives 



SVix) 

d_ 

dx 



(x) {g{x)t{x) - 
--^Pv{x)v'{y)g{x) 



[P'(^).9'(x)]' 



(20) 



According to the statement above, letting g{x) = 1, we 
instantly have 



1 



t'{x) = -7;Pv{x)v'{x) 



Integrating on both sides, then 

1 poo 



t{x)dx 



(21) 



dx I dx' py{x^)v' {x^) 

) >/ —OQ 

/•OC 

dx' I dxp^{x')v{x') (22) 

) J x' 

x' py{x')v'{x')dx' , 



which coincides with the virial theorem 6]. Stronger 
than the virial theorem, Ea.l|21l) provides another 
functional relation (see e.g. Baltin Jj, March and Young 
[3 to be satisfied by the pair of kinetic energy density 
and density so that they can possibly be v-representable. 



III. PRIMITIVE APPLICATIONS 

Due to the last nonlinear term on the left-hand-side in 
the constraint Ea. (|17|l . it can not easily be simplified. As 
to its positivity, by making the constraint a bit stronger, 
we can eliminate this term so that the feasibility of the 



constraint can be illuminated(Empirically, the correction 
due to the {g')^ term always turns out to be very small 
- for an exception, see Ea. (|32|l '). Minimizing the left- 
hand-side with respect to the transform g{x), we have 
the constraint on t(x) and p{x) simplified as 



Rv[p,t] 






1/3 



dx- 



pV{f)dx>EM[V], (23) 



where ^p and / are functions of x, and 



V 



'pV'if) 






1/3 



(24) 



The ordinary differential equation(ODE) array above 
won't challenge numerical calculations at all. However, 
in order to present the constraint Ea. H23|) concretely, let's 
consider a reference system with a linear external poten- 
tial V{x) in the one dimensional half space. That is 



V{X) = 



X if a: > 0, 

oo otherwise 



With this, the ODE array fEa l^ is instantly solved, and 
the constraint can be further simplified to 

Rl[pA^IJ {A{x)f/\2t{x)f/^dx>EL{N), (25) 

where A{x) is the cumulative density, defined as 



A{x) 



p{x')dx' , 



(26) 



El{N) ~ ^ Si=i 2^^^'^ai, the ground state energy for 
N non-interacting fermions with a linear well as exter- 
nal potential, Oj is the ith root of the first kind of Airy 
function Q- As the number of fermions N increases, 
asymptotically, we have El{N) « ^{i-nf/'^N^/'^, with 
extremely rapid convergence. 

If we are interested in the asymptotic value of the 
ground state energy, we can apply this to finding the 
ground state energy of a non-interacting fermion system 
with half space external potential v{x) — —x^' 



mini?[pjv] 

Pn 



. 7 + 2 
mm 



PN{x)x'^dx (27) 



with the constraint 



Pjv 27 

RL[p,t] = \j A^'^{2tY'^dx > l(37r)V3ivV3, (28) 

where the local virial constraint Ea. (|21|l has also been 
applied. It can be shown that the minimum of Ea. H27|l 
occurs only at the boundary of the function space of the 
density pn{x). And since the density function is non- 
negative, we have pn{x) = NS{x — xq). From Eq.(|2Hl, 
we have 



xo>{l 



3 

y+2 



{SttN)-- 



(29) 



Therefore, 



Substituting a;o from Eq|21| then 



E[pn]/N > 



7 + 2 
27 



3-1 

T + 2 



(SirN)^ 



Not surprisingly, when 7 = 1, we have 

E[pn]/N > ^(3^7V)2/3, 



(30) 



(31) 



which is identical with the exact asymptotic behavior of 
the system with the half space linear well as external 
potential. For the Coulomb potential, where 7 = — 1, 
Ea. (|Sn|l asserts that 



125 

Ec[pn]/N>-j^N- 



(32) 



which is clearly false because Ec'[pn] — > - const as 
N —> 00. The reason for this apparent paradox lies in 
the neglected (g')^ term, which is ordinarily very small, 
but when it is required to map a density due to a reg- 
ular potential onto one from a singular potential, this 
is no longer the case, and the inclusion of (g')"^ could 
be mandatory. We have taken the 7 = — 1 case as a 
first example to show that the reference system must be 
reasonably similar to the system under study to make 
sense. Ignoring the (g')^ contribution accentuates this 
difference to the point that a lower bound is no longer 
obtained, but the general comment remains solid. We 
now consider further examples that develop this implicit 
criterion. 

For harmonic oscillating fcrmions with 7 = 2, we have 



Eho[pn]/N> (- 



3/2 



3t:N « 0.842978iV 



(33) 



compared to the exact l.OiV. 

As 7 ^ 00, we have the rigid wall box to solve, the 
constraint gives 



Ensip.] > ^^N^ 



(34) 



6 



compared to the exact asymptotic Ej.^„^^ u _^ 
As we see, with one constraint from the linear well ref- 
erence system alone, applying the minimization scheme 
won't result in a bound very close to the exact solution. 
If we subject the density minimizing the energy func- 
tional above to a different constraint, for example taking 
harmonic oscillating fermions as a reference system, we 
have 



Rho [p, t] = min / 

fix) \Jo 



N 



2in^)y- 



;dxXip ^ 



--Nx, 



N 

1+2 



ifi^o))'- 



(35) 



RHo[p,t]^N 



3/2 



(37r7V) « 0.847V2 (36) 



Thus, for the linear constraint optimum, constraint 
RHo[p,t] > N^ is badly unsatisfied. As will be shown 
later, the ground energy level obtained with the min- 
imizing scheme will be dramatically improved if more 
reference systems are involved. 



IV. DENSITY CONFIGURATION 
CONCENTRATION 

For any A^-representable density profile p{x) , there al- 
ways exist scale transformed density functions pa{x) = 
ap{ax) so that using virial constraint (|22|l the applied 
reference dependent constraint is satisfied. Among these 
transformed densities, we have a uniquely defined refer- 
ence dependent density functional 



E^ 



Rem = ™n 

a satisfying 
the constraint 



KEDjjY(a;; [pa]) dx 



+ / Pa{x)v{x)dx 



(37) 



where KEDjjy is the kinetic energy density term with 
the local hyper-virial constraint H21|l applied. Therefore, 
a pseudo energy landscape (PEL) of the density profile is 
obtained, the global minimum of which provides a lower 
bound to the energy expectation of the ground state. It is 
noticed that because of making use of the scale transform, 
the landscape -£'j^gf[p] depends only on the scale free 
configuration of the density profile. Rather than search- 
ing the minimum over the entire non-negative density 
function space, one can obtain the identical minimum by 
searching within a more compact subspace provided that 
every configuration of density has been included. There- 
fore, the reference dependent landscape is furthermore 
reduced to the density configuration hypersurface. The 
minimization task can be achieved by searching the mini- 
mum over the randomly generated n dimensional discrete 
density profile subspace (100000 configurations were typ- 
ically used ) as 

Po ® Pi "Xi • • • (Xi Pn-i ^ 0® Rand{m) (8) • • • ® Rand{m) ®Q 

^ V ' 

«-2 

(38) 
where both n and m are positive integers and Rand{m) 
randomly generates an integer between to m — 1. po and 
Pn-i are zero since the density vanishes at both ends in 
one dimensional space. The entire configuration hyper- 
surface is covered as n and m approach infinity. To avoid 
the unnecessary computing effort in searching within the 
unsmooth function space, a smooth density function sub- 
space can be obtained by interpolating the randomly gen- 



erated control points with splines. Therefore, for prac- 
tical use, finite number of random numbers is required 
to generate a density profile. In this paper, for random 
generated density profiles, a B-spHne [i3 of degree 3 is 
applied. 

As shown in Figure [21 the control points can be ei- 
ther equally spaced or randomly spaced to include more 
rapidly varying density configurations with the same 
number of control points. This is useful when, as in this 
introductory study, we want to economize on the resolu- 
tion employed. 

The lower bound of the ground state energy is given 
by 

El ^ min ^ppf [p = spline (pi, p2, •■ -Pn-i)] 

Pl,P2,■■■Pr^-l 

(39) 
According to the definition of the landscape -Ejjpf, if 
one were to take reference system with external poten- 
tial identical to that of the target system, it's trivial that 
the reference dependent landscape -£'j^gf[p] = -Egxact' 
where inexact i^ ^he ground state energy of the target 
system, which is flat all over the density configuration 
hypersurface. Therefore, any landscapes tangential to it 
have the landscape altitudes concentrating near the ex- 
act ground state energy, which is a saddle point. Note 
that for a reference system very different from the tar- 
get system, there's only one crosspoint in the density 
configuration hypersurface between the reference depen- 
dent landscape and the fiat landscape described above. 
It can be concluded that the distribution function of 
^Ref ^^^ ^^^ maximum at the ground state energy of 
the target system, provided none of the density config- 
urations dominate when generating the density profiles 
randomly, as illuminated by Figure |3| Consequently, the 
randomly generated density configurations provides not 
only a lower bound but also practical approximation to 
the target system. 



TABLE I: Ground state energy of Af-harmonic-oscillating- 
fermion obtained from reference system with linear well. 
-Eoxact is the exact energy expectation solved for N non- 
interacting fermions, and -Emost is the energy expectation with 
highest concentration generated by random density profiles 



N 



Sc: 



E„ 



1.5 

5 

10.5 

18 



1.6 
5.0 
10.5 
17.8 



exact energy expectation 
'most probable energy 



Consider the non-interacting iV-fcrmion system with 
v{x) — ^ as the external potential to solve. Rather than 
taking the asymptotic constraint, one can make use of the 
exact energy level of N non-interacting fermions within 
the linear well. The pseudo energy landscape is generated 
with the -Epjgf associate with the random density profiles. 



which are B-spline interpolated. We find that the energy 
distribution concentrates at the following values shown 
by Tabled The exact value is provided as a comparison. 



V. MOLECULE-LIKE SYSTEMS 

In our primitive application of the energy minimizing 
scheme, we have shown that coordinate-transformations 
from a known reference system do provide a somewhat 
weak lower bound for the ground state energy level of 
a fairly different system, too. We are also interested in 
applying the constraint to a locally similar but globally 
different target system . Let's see what would happen if 
we convert a single well reference system into a molecule 
with double core. 

For example, let's consider the fermion system to be 
solved as having the double linear well as external poten- 
tial 



v{x) 



X -~ b ii X > b, 
b- X iiO < X <b 
x + b ii -b < X < 
, —b — X if a; < —6 



(40) 



The reference system has the external potential illus- 
trated by Figure 21 



Vix) 



a{x ~ b) if X > b, 
a{b ^ x) if X < b 



(41) 



By minimizing the energy subject to the single well 
constraint, the minimization take place at the boundary 
of the function space of density. 

N 



Pn{x) = y (<5(a; 



xo) + 6{x + xo)) , (42) 

and Xo is determined by the asymptotic constraint 



Rls[pn] ^\ j {A,{x)f'\2T{x)f'^dx 



> 



10 



2/3 



(43) 



Af5/3 



, where 



Asix) 



J^ PN{x)dx ifx>6, 
/_ pN{x)dx ii X < b 



We have 



a:o>-l 



1 /37riV\ 

2~y 



2/3 



(44) 



(45) 



And so we have for the ground state energy of the 
double linear well system 



E[pn] = -^Nxo 



Nb> N 



3 /'3TrN\ 

iov^~; 



2/3 



(46) 
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FIG. 2; (Color online)To illustrate how the random density configuration is generated, the solid lines show the random density 
profiles, which are provided by B-splines along the arbitrary selected set of randomly generated control points marked with 
'o's. (a) shows the instance with 11 control points evenly spaced while (b) shows the instance with 11 control points randomly 
spaced. The density configuration space generated with evenly spaced control points is a subset of density configuration space 
generated with the same number of randomly spaced ones. 



as a lower bound. 

On observing the energy landscape concentration as 
shown before, we find the most probable energy levels, 
which can also be compared with the exact ground state 
energy level of the double linear well system. As shown in 
Tableim the most probable energy levels give the best ap- 
proximation. With N increasing, the minimization with 
single well constraint alone gives a closer result than the 
Thomas-Fermi approximation. 

The potnetial consequence of the quite decent agree- 
ment at this preliminary level are far - reaching: the 
same strategy can be used for full interacting atoms in 
molecules, a situation in which universally useful model 
systems are few in number. 



VI. MORE THAN SINGLE CONSTRAINT AND 
SHIELDED COULOMB POTENTIAL 

As stated above, since the density configuration for the 
system to be solved should agree with every reference sys- 
tem provided, the more constraints we apply, the better 
the result one would expect. Therefore, we would like to 
constrain the target system by two reference systems to 
see how the lower energy bound is improved. 

Among all the density functions which satisfy the con- 
straints, we are interested in those that generate the low- 
est energy levels together with the assistance of the local 
hyper-virial theorem. These energy levels are used to 



TABLE IL Ground state energy of double linear well system 
with 6=1 and single linear well reference. The second column 
is the exact energy for non-interacting fermions with double 
well as external potential. The third column is the energy at 
which the density configurations concentrate with a single well 
as reference system. The fourth column is the energy obtained 
by subjecting the minimization to the single well constraint, 
serving as lower bound. The fifth column provides the result 
given by Thomas-Fermi Approximation as a comparison 



N 


-Ecxact " 


-Bmost single well' 


-Emin single well " 


-Bt-f 


1 


0.6266 


0.630 


-0.157 


0.5312 


2 


1.6622 


1.580 


0.677 


1.6118 


3 


3.5386 


3.206 


2.262 


3.1791 


4 


5.9278 


5.477 


4.499 


5.2861 


5 


8.9200 


8.337 


7.328 


7.9197 


6 


12.421 


11.74 


10.71 


11.059 


7 


16.430 


15.67 


14.60 


14.684 


8 


20.916 


20.08 


18.98 


18.775 


9 


25.861 


24.95 


23.84 


23.317 


10 


31.252 


30.28 


29.14 


28.295 


11 


37.069 


36.05 


34.88 


33.696 


12 


43.305 


42.23 


41.04 


39.508 


13 


49.944 


48.82 


47.61 


45.720 


14 


56.979 


55.81 


54.57 


52.323 


15 


64.399 


63.19 


61.93 


59.309 


16 


72.196 


70.94 


69.67 


66.668 



"exact energy expectation 

''most probable energy with single- well constraint 

"^lower bound with single-well constraint 

''energy expectation with Thomas-Fermi approximation 



provide a lower energy bound. The B-splines generated 
by random numbers are still used as the density configu- 
ration candidates. 

We will see that not all pairs of constraints are effective 
in providing a decent lower energy bound for the target 
system. Not surprisingly, only those pairs of constraints 
that bracket the target system are really effective. 

Note that only the constraints from the reference sys- 
tem with external potential v{x) ^ x can be written down 
explicitly, while the others we can make use of so far can 
only be applied numerically. A primitive test of the effect 
of the double constraint is by making use of the asymp- 
totic behavior of the constraints when the number of par- 
ticles is sufficiently large, and comparing the result with 
those obtained from Thomas-Fermi approximation. The 
latter is believed to be a good approximation to describe 
the energy levels for the non-interacting fcrmion system 
when the number of fermion is large enough, except for 
those with Coulomb potential as external potential. 

The target system we will test has external potential in 
the form v{x) = ^. Pairs of constraints that can bracket 
the target system and those that can not are selected to 
generate the lower energy bounds for the target system. 
The results are shown in Table ITTll We see that the re- 
sult doesn't improve much compared with those obtained 
from the single constraint unless the pair of references 
system can bracket this system. The role of "bracketing" 
is overwhelmingly important. 

As discussed above, we have shown that our continu- 
ous coordinate transform scheme works well with most 




"exact 



£ref[P] 



^ 



Ref, 




P 



FIG. 3: (Color onIine)The cartoon shows the landscapes E-n^^ 
in density configuration space, which have dependence on dif- 
ferent reference systems. For those reference systems with 
similar external potential to that of the target systems, such 
as Refo shown in the figure, we have a fiat landscape rep- 
resented with a horizontal line. For reference systems very 
different from the target system being applied, the reference 
dependent landscapes are shown in the as dotted curves. All 
the landscapes concentrate at the exact ground state energy 
level of the target system. 
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FIG. 4: (Color online)An illustration of the potential of ref- 
erence system and target system. 



TABLE III: Lower energy bound obtained with constraints 
from two reference systems. The asterisk in the first column 
indicates that the pair of reference systems can bracket the 
target system. 



7 


Ref'" 


Ref 




-C/most 


-Et-f 


2' 


X 


x' 


0.99A/-^ 


1.03Af^ 


l.OOA/-" 


2* 


-x-^ 


x' 


0.64iV2 


1.02iV2 


l.OOAf^ 


2 


1 


a;3 


0.37iV2 


/ 


l.OOAf^ 


2 


-x-^ 


X 


0.84iV2 


l.OOiV^ 


l.OOiV^ 


2 


1 


X 


0.84iV2 


/ 


1.00A/-2 


1.5* 


X 


x^ 


1.09Ar^ 


l.lliV^ 


1.107V¥ 


1.5 


X 


1 


l.OOA/-^ 


/ 


l.lOAf^ 


3 


X 


2 
X 


0.83iV^ 


0.93Af^ 


0.94Ar^ 


3 


X 


1 


0.68A/-Tr 


/ 


0.94Ar^ 



"the external potential of the applied reference system 
''lower energy bound with the constraints applied 
'^niost probable energy with double constraints 
''asymptotic energy expectation with Thomas-Fermi Approxima- 
tion 



fermion systems having an infinite number of bound 
states. It will be interesting to know if it will break 
down when the system has finitely many bound states. 
One such A^-fermion system that was the focus of past 
study is the shielded Coulomb system 11] with the exter- 



nal potential v{x) 



r/D 



, which can only be numer- 



ically solved d2||. It's known to have a finite number of 
bound energy levels due to the effect of screening, and the 
number of bound energy levels depends on the screening 
length D. 

The pair of reference systems we use here are Coulomb 
and harmonic oscillator systems. And the exact value of 
each energy level is applied for each reference system. 
For half space 1-D Coulomb systems, the ground state 
energy level for A^ fermions is proportional to H{N, 2) = 

12n=i ~^' where H{x, 2) is the harmonic number function 
of X with order 2. And the ground state energy level 
for A^ fermions within a half space 1-D harmonic well is 
proportional to A^(A^ +5)- 



The results with random density configuration candi- 
dates are shown in Table IIVI We do have finite num- 
ber of bound states with this pair of constraints, even 
though the reference systems that we apply have an in- 
finite number of bound states. Only when the screening 
length arrives at some threshold length, can a new bound 
state survive from the double constraints. 



TABLE IV: Lower bound for the ground state energy with 
shielded Coulomb potential. The number to the left of each 
cell is the lower bound obtained with double constraints. The 
number to the right is the energy level obtained with numer- 
ical calculation. 



D 


1 


N 
2 3 


4 5 


2 

5 
10 
20 


-0.244|0.148 
-0.397|0.327 
-0.4480.407 
-0.4710.452 


/ / 
-0.41910.339 / 
-0.521|0.457 -0.52610.460 
-0.571 0.534 -0.602 0.553 


/ / 

/ / 

/ / 

-0.607|0.556 / 



VII. CONCLUSION 

We have seen that the analysis of a many-fermion 
ground state can be recast as a constrained minimiza- 
tion of a functional of two scalar and one tensor field. For 
non-interacting fermions, only the density and a suitably 



defined kinetic energy tensor density are required, and 
the effort is shifted to tabulating and using the needed 
constraints, which can be tuned to emphasize known and 
hypothesized physical aspects of the system. Taking the 
one-dimensional spinless system as prototype, a class of 
inequalities based upon solvable models has been devel- 
oped, as well as the hypervirial equality. With only a 
small number of solvable models to help us, quite de- 
cent results have been obtained for some simple systems, 
including mock "molecular" ones and systems with a con- 
trolled number of bound states. Extension of these tech- 
niques to 3 dimensions, and to the spin degree of free- 
dom, is direct, and will be reported in the near future. 
Extension to physical interaction is less direct, and sev- 
eral paths are under study, which will be reported as well. 
They include universal bounds on the interaction energy 
as a functional of p and t (see refQ for a rudimentary ex- 
ample), and modification of the coordinate transforma- 
tion tactic in the face of interactions. It must of course 
be emphasized that this continuing investigation is not to 
be regarded as competition to the array of high accuracy 
computational techniques that have been developed, but 
rather as a low-cost replacement when fine detail is not 
required. 
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